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Abstract 

We investigate the sigma and baby Skyrme models with an RP^ target 
space. We compare these models to models with an target space. We 
investigate the interactions between solitons and defects in the RP^ sigma 
model. 



1 Introduction 

The Skyrme model has long been of interest as an effective field theory of nu- 
cleons [0. The 2 dimensional baby Skyrme model allows us to study a more 
tractable analogue of the Skyrme model §. The related 0(3) sigma model is 
applicable to some condensed matter systems. Nematic liquid crystals are a 
condensed matter system where the target space is not but RP^. In this 
paper we describe our work on the variants of the baby Skyrme model and 0(3) 
sigma model using MP^ as a target space. 



1.1 Models 

The 0(3) sigma model describes maps from (the "physical" space) to S'^ 
(the "target" space). is the normal Euclidean plane, whereas S"^ is the 2 
dimensional surface of a unit sphere embedded in 3 Euclidean dimensions. The 
lagragian of the 0(3) sigma model is given by 

C^^df^^-dJ, (1) 
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where cf) gM.^ such that (p-cfi = 1. The condition that 4> '-^ 4'vac as r i-^ oo, where 
r is spatial radius, is usually imposed to keep the action finite. This effectively 
compactifies the physical space from to S^. This in turn implies that the 
configurations must fall into disjoint classes characterised by the elements 
of the second homotopy group 7r2(S'^) = Z. Field configurations which fall 
into any of the classes corresponding to non-trivial elements of 7T2{S'^) describe 
'lumps' of energy. These lumps are not stable, and as a result are not true 
solitons. This instability is due to the conformal invariance of the model - the 
lumps have no intrinsic scale, and so may change scale without any energetic 
penalty. This scale invariance may be removed by the addition of new terms to 
the Lagrangian, as seen below. 

The first model we examined was the baby Skyrme model with the so called 
'new' potential term. The Lagrangian is 

£ ^ ^d''^ ■ dj + 0i^{{d>'$ ■ d''$){dj ■ dj) ^ {d''^ ■ djf) + e2Cv, (2) 

where 9i and 62 are real positive constants. The term containing 61 is the 
so called (2 dimensional) Skyrme term. This term tends to cause solitons to 
broaden, and is the only possible Lorentz covariant term of fourth order with 
only first order time derivatives. The term containing 62 is normally refered 
to as the potential term. This term is chosen so that it acts to shrink the size 
of the soliton. Consequently a stable equilibrium may be reached where the 
effects of the Skyrme term and potential term balance. A number of different 
potential terms have been studied ^ , but our MP^ target space forces us to 
use a potential term which is invariant under the transformation (j) i-^ —<j). A 
good choice would be the so called 'new' potential 

Cv^\{l-{$-$.acf). (3) 

1.2 The Geometry and Topology of MP^ 
1.2.1 Geometry 

RP^ may be thought of as a hemisphere, and so may be parameterised by a 
unit 3- vector (p where = — RP^ may also be described locally by a 
2-vector x where 

x; = ^ Z-1..3 j = 1..3. 

Here i labels which chart is being used, and j labels the components of x- X 
appears to have 3 components in this definition but, as the component x\ is 
always 1, x has only 2 degrees of freedom. The chart labeled i is ill defined on 
the equator = 0, and so all three charts are required to cover RP^ - the first 
covers all of the manifold except for an equator, the second all of this equator 
except for 2 antipodal points and the third those two points. 
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A neat way to combine all three x charts is to use the matrix Pij [x) = 
4>i{x)(j)j{x). Each chart is easily extracted as = -p^ with no sum implied by 
repeated indices. 

1.2.2 Topology 

The nth homotopy group of KP^ can be found quickly for n > 2 as 5"^ is a 
universal covering space of RP^ . A theorem found in Q states that this implies 
that 7r„(S'2) 9^ 7r„(RF2) V n>2. 7ri(Rp2) may be found to be the second 
symmetric group, Z2. To summarise: 

'Ki{S'')^e 7ri(Rp2)^Z2, (4) 
7r2(S'2)=Z 7r2(Rp2)=Z, (5) 
7r3(S'2)9^Z TT-i{m.P^)^'L, (6) 

1.2.3 Winding Number 

The class of field configuration may be characterised by the degree of map, 
or winding number. To calculate the winding number we must calculate the 
puUback of the two form on the target space; the map from the physical space 
to the target space induces a natural map (the "puUback") from the two form 
on the target space to a two form on the physical space. If the two form on the 
physical space is then integrated over all space the result must be a multiple of 
an integer. By renormalising we get a formula for the integer valued winding 
number: 

J cfx€ijeabc4'aidi(j)b){dj4'c). (7) 

It should be noted that the transformation cj) i— > —0 leads to T — T. This is 
because RP^ is a non-orientable manifold - a two form is not well defined on 
RP^ as its sense is changed by a translation along a non-trivial loop. In fact it 
may be shown (s) that field configurations of (f are characterised not by 7r2 {A4 ) 
but by 7r2 ( ) /tti ( ) . In the case of an RP^ target space this means that field 
configurations are characterised by Z/Z2. 

1.3 Numerical Techniques 

The models outlined above are highly nonlinear and more or less completely 
intractable analytically. The models are, however, open to study using numeri- 
cal techniques. One may write a simulation on a discrete grid and analyse the 
evolution of an arbitrary initial state. With a clever ansatz for the initial condi- 
tions, with or without the use of a relaxation routine, the behaviour of minimal 
enery field configurations may be studied. 

We simulated the sigma and Skyrme models numerically, using RP^ as the 
target space. Our simulations were based on a 200 x 200 point grid using the 
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nine point laplacian and derivatives outlined in Q , together with a fourth order 
Runge-Kutta algorithm for simulating time evolution. The timestep length was 
set to half the gridpoint spacing. To give the M.P^ topology the state vector was 



stored using the P matrix described in section (1.2.1) for position on the target 
space, and Qij = Pij = (j)i4'j + 4'i4'j for the rate of change of position on the 
target space with time, where a dot is used to denote derivatives with respect to 
time. The simulation then evaluated derivatives by determining the field as the 
equivalent (f> for any given set of nine points, with the central point mapped to 
the north pole of S"^ . Reflections from the boundaries of the grid were reduced 
by damping out kinetic energy in a region near the boundary. The fields were 
kept on manifold by applying the following transformations each timestep: 

^kk 



These transformations ensure that Tr{P) = 1 (i.e. the field lies on MP^) and 
that Tr{Q) — (i.e. the rate of change of the field is tangential to the surface 
of the field manifold) . 



2 The Valued New Baby Skyrme Model 

Simulations of this model were in close qualitative and quantitative agreement 
with simulations using an target space. 



2.1 The Hedgehog Anzatz 

As in the S'^ model the simplest ansatz for a skyrmion is the radially symmet- 
ric ansatz, usually refered to as the 'hedgehog' ansatz [Q. For a skyrmion of 
topological charge one this ansatz gives: 

/ sin[/(r)] cos{9 - 7) \ 
$^ sin[/(r)] sin((? - 7) . (10) 
V cos[/(r)] ; 

Here f(r) is known as the profile function and 7 is an arbitrary phase parameter. 
This profile function is arbitrary up to the boundary condition that /(O) = 
(2n + l)7r and /(oo) = 2m7r where n,m € Z. To minimise the energy of the 
ansatz one may determine the energy of the field as a functional of the profile 
function and its first derivative, as shown below: 

/•oo 2 r 

E = 2n rdr{f'^ + ^(1 + 26^/^) + e^V{f)). (11) 

One may then use the calculus of variations to extremise this energy with the 
appropriate limits at r = and r = 00. The resulting Euler-Lagrange equa- 
tion is a non linear second order ODE for the profile function. We solved this 
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numerically, using the shooting method. This ansatz leads to a stable single 
skyrmion. Note that our ansatz is identical to that of the skyrmion, and so 
our result agrees with the result of Weidig Q . Simulations showed this ansatz 
to be stable, again in agreement with Weidig. Note that the energy has been 
defined such that the energy of a hedgehog ansatz soliton with 6*1 = 02 = is 
given hy E = 1. 



2.2 Two Skyrmions 

As for the model, a two skyrmion ansatz may be arrived at by taking the 
stereographic projection from the north pole of a one skyrmion: 

= ^ 7 . (12) 



This may then be combined with another skyrmion using the ansatz 

1 _ 1 1 



(13) 



so that the final field Wt is approximately equal to the field for any skyrmion 
near that skyrmion (where W = 0) and the final field takes the vacuum value 
of VF = oo at spatial infinity. Using this ansatz we were able to reproduce 
the attractive and repulsive channels found in Q, as well as the ninety degree 
scattering found in the attractive channel. 



2.3 Equivalence of S"^ and RF^ models for Smooth Maps 

In fact the two models are exactly equivalent as may be seen from the following 
argument. Consider maps /, (j): 

f -.r^ 5", (14) 

where n < m, and 

: S"" M. (15) 

As 7r,i(S'™) = e V n < m f may always be smoothly deformed to map all 
of /" to a point. This implies that, for all smooth maps <j), all ?i-loops in 5'™ 
shall be the pre-image of an n-loop in M corresponding to the identity element 
of 7r„(Af). 

One consequence of this is that for M = MP^, n — 1 and m = 2, all maps cj) 
may be described as maps to M = 5^, as RP^ is universally covered by S^, so 
the RP^ new baby Skyrme model is identical to the 5^ baby Skyrme model if 
(j) is well defined at all points. 

It is also worthy of note that this argument suggests that simulations of 
Fadeev-Hopf solitons [Q, |j using RP'^ as a target space should reproduce the 
simulation of such solitons with an target space unless defects are present. 
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Figure 1: Maps from /„ to Sm to M 

3 The Sigma Model with Defects 

The implication of the argument put forward above is that to discover new 
behaviour in an MP^ model we must introduce a discontinuity into the field. 
For such a point to be stable a circle around the discontinuity in physical space 
must map to a non-trivial loop in the target space. If this is not the case then 
the discontinuity can be removed by a continuous variation in the field. These 
(point) discontinuities are analogous to disclination lines in a liquid crystal. We 
shall refer to them as "defects" . One consequence of having an odd number of 
defects in the system is that a contour around the edge of the physical space 
must map to a non-trivial loop on EP^ . This means that it is no longer possible 
to have the field tend to a single value at spatial infinity. As a result, it is no 
longer easy to use a conventional potential term such as (^) in the Lagrangian. 
We may use the KP^ sigma model, so our model has neither Skyrme nor a 
potential terms, or we may use an unconventional potential term. One possible 
potential term would be of the form 

Cv = ^(0 • $mass)'^), (16) 

which we shall call the "easy plane" potential. This effectively gives a mass 
to one of the fields {(j)mass) where the potential (||) gives mass to both fields 
orthogonal to (pvac- Easy plane baby skyrmions have not previously been studied 
and so are worthy of attention in their own right. We have confined our work 
to the sigma model in all that follows. 

It should be noted that if there is no single value for the field at infinity we are 
no longer compactifying the physical space to 5*^. As a result winding number 
is no longer necessarily conserved. In fact, the physical space is topologically 
iS''^ X M} as not only is the space not compactified at infinity, it also has no well 
defined map at the point of the defect. Whilst the mapping from the target 
space is frozen at infinity, it is not frozen around the defect. The circle around 
the defect necessarily maps to a non-trivial curve in MP^. If we represent KP^ 
as a sphere, where antipodal identification is allowed, a non-trivial curve is a 
line between a pair of (arbitrary) poles. If we have a defect-soliton system, one 
of these curves represents both the image of a circle at infinity in the physical 
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space, and the image of an infinitesimal circle around the defect. As a circle 
is contracted in from infinity to the defect, the image curve rotates around the 
sphere, hinged on the polar identification, wrapping the sphere once. As the 
image of the infinitesimal circle around the defect is able to move, it may rotate 
around the identified poles, unwinding the soliton like object. 

We decided to study the behaviour of a soliton like lump in the presence of 
a point defect. To do this we need an ansatz for a soliton-dcfcct system to use 
as an initial conditon for our simulation. Such an ansatz needs to be sufficiently 
close to equilibrium to avoid large quantities of radiation perturbing the system 
- if the system is far from equilibrium much of the excess energy will be shed 
from the solitonic objects as radiation. 



3.1 Defect Ansatz 

The most obvious ansatz for a defect is a radially symmetric, planar field. A 
defect at r = may be written in polar coordinates as 



, (17) 

:osf 

which, using the stereographic projection from the south pole, 

W=^±^, (18) 



is given by: 



W" = tan-. (19) 
4 



This ansatz is defined for < ^ < 2-k. On the line ^ = there appears, at 
a first glance, to be a discontinuity in the field, but as 4> is identified with —0 
there is no such discontinuity. In fact any field configuration which contains a 
defect must have a line from the defect to infinity along which the identification 
between and —(h is made. It is simple to show analytically that the above field 
configuration is in equilibrium. Simulations also show that this field configu- 
ration is stable. It should be noted that this field configuration has an infinite 
action arising from a logarithmic; divergence as r tends to infinity or as r tends 
to zero. This is not a problem if the physical system we are modeling is a liquid 
crystal. This is because such a system is, firstly, finite in extent, and, secondly, 
a discrete system. The field approximation is clearly not appropriate at r = 0. 
Care must be taken when simulating such an object numerically to place such 
an object away from a grid point. 



3.2 'Glued' Ansatze 



7 



3.2.1 The Three Region Ansatz 



One ansatz for a defect-soliton system which we tried early on involved 'glueing' 
together three regions in an attempt to graft together a known ansatz for the 
soliton to a defect like object - see figure (0). In region I the field followed the 



Region I 



Region n Region III 



Figure 2: Regions and Coordinates for Glued Ansatze 

standard lump ansatz for the 0(3) sigma model. This is given by the hedgehog 
ansatz (^|) with the profile function given by: 

9\r 

/M=tan-iY3AV^, (20) 

where A is an arbitrary (real) parameter describing the (inverse) 'width' of the 
lump. In region II the field is constant and equal to the vacuum value of region 
I - there should be little discrepancy at the border between regions I and II as 
a result, and any variation may be interpolated across several grid points in the 
simulation. In the region III the field is given by: 

/ cosOd \ 

$={ , (21) 
V sin^rf J 

for —'!T<6d<Tr. The region II and region III fields match along the line 
= ^ and, due to the identification (/>=—</), along the line 6 = This 
places a defect at the origin. Unfortunately this defect is far from stable - the 
energy of this defect is double that of the defect described by dlj). Thus in the 
simulations of this total field configuration we have seen the emission of waves 
of radiation by the defect. These waves affected the skyrmion. To reduce these 
problems we have modified our ansatz as detailed below. 
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3.2.2 The Two Region Ansatz 

We modified the above ansatz by replacing the field in regions II and III by the 
field below: 

W = tan ^e-'=(''''-''«\ (22) 

where < dd < 2tt and W is the stereographic projection from the south pole. 
At a first glance this field looks as though it might tend to the vacuum value 
in all directions whilst looking like a defect at small radii. In fact, this is not 
the case - the field is smooth away from the origin, and therefore a loop around 
this field at infinity must be non-trivial. The field must vary dramatically as 9 
approaches 27r and r approaches infinity. We had hoped that this would not be 
a problem in the (finite) region of the simulation, but discretisation brought its 
own problems, as shown in figure (||). The line where W — I must cross the last 
line of grid points on the large dd side of the dd = line. This leads to a point 



w~o 




Figure 3: Diagram showing discretisation problem 

where the variation in the field from one grid point to the next is as large as it 
can be, resulting in a region of high energy density. In fact this region also has 
defect type winding number - the ansatz has created a second defect in a far 
from equilibrium configuration. In fact if one constructs an ansatz which tends 
to constant field at spatial infinity (or the edge of the simulation region) then 
the space must contain an even number of defects. A loop at large r which maps 
to one point in M.P'^ must correspond to the trivial element of 7r2(MP^) = Z2, 
so the number of defects inside this loop must be even. 
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3.3 The Stereographic Defect - Soliton Ansatz 

To create an ansatz for a field configuration including both a soliton and a defect 
we must re-express our field configurations using the stereographic map, W. A 
projection from the south pole causes the north pole, (/)3 = 1, to map to = 
and the south pole to map to ly = cx3. Note that the antipodal identification 
under this map is given by 

W^--^. (23) 

W* 

Our ansatz must a) look like a soliton near the soliton, b) look like a defect near 
the defect, c) look like a defect at infinity and d) be smooth at all points away 
from the defect. 

A single soliton at the origin may be described by the field 

a cos{9 — 7) 

asin(6l-7) | , (24) 

where a — j^jrpi, P — , 7 is an arbitrary (real) phase parameter, A is 

an arbitrary (real) parameter characterising the width of the soliton, and (r, 9) 
are polar coordinates. Note that this is the hedgehog ansatz with the profile 
function given by (|2^). To combine this with a defect we choose the point of 
projection of both fields such that the soliton at infinity maps to a number of 
modulus 1, and the line along which the defect field identifies approaches from 
one side and 00 from the other. These fields are given below. 

^ (25) 

1 + acos(t' — 7) 



1 + cos(|) 





tan(-). (26) 



If we multiply these fields we find 



T^. = ^^^^fc^tan(^), (27) 
1 + q;cos(6's - 7) M ' ^ ^ 

where 9d is the polar angle coordinate with respect to the position of the defect 
and (r, 9s) are the polar coordinates with respect to the position of the soliton. 
This ansatz obeys conditions c) and d) above, c) is satisfied as at infinity 
Wt = —i X Wd- This is acceptable as multiplying a field configuration by a 
number of unit modulus is merely a rotation about the axis of projection, d) is 
also satisfied as the line of identification in Wd is preserved by the multiplication 
-Ox Ws ~ and (crudely speaking) 00 x Wg = 00. Condition a) is fulfilled 
if the soliton and defect are sufficiently separated for the variation in Wd to be 
small across the scale of the soliton. Condition b), however, is only fulfuUed if 
{9s - 7) = (2n + l)f where n € Z. 
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3.4 The Inhomogeneous Defect - Soliton Ansatz 

An alternative ansatz involves using the inhomogeneous coordinates outlined in 



section (1.2). If we define a complex number 



(pi + i(l)2 

W = ^ , ^ (28) 

03 

then we can express the soliton field as a complex number which tends to at r 
tends to and oo. If we add this to the defect field then the result must satisfy 



conditions c) and d) outlined in section (3.3). c) is satisfied because the map 
used is a map of MP^, not 5^ and so needs no line of identification for the defect 
ansatz, whilst d) is obeyed as the soliton field vanishes at infinity, leaving only 
the defect field. Explicitly, the fields look like this: 

m = ^e^(^-^), (29) 
P 

Wd = tan(^), (30) 

WT = ^e^(^=--)+tan(|). (31) 

This expression obeys a) if Wd is small around the position of the soliton. 
Condition b) is only met if Ws is small in the region of the defect, and in 
the example above condition b) is less well satisfied if Ws has an imaginary 
component at the defect. 

3.5 Simulation methods 

For simulations of a defect-soliton system we introduced a conformal grid in 
a similar manner to Leese et al in |9j. This involved changing the physical 
coordinates {x,y) to {x',y') such that 

This allows a much larger area to be simulated. The grid was 199 x 199 points 
with dX' = 0.01. This conformal grid was introduced to reduce boundary 
effects. Again, the timestep length was set to half the gridpoint spacing. 

Another alteration to the simulation for defect-soliton systems was the in- 
troduction of a period of relaxation before the simulation was allowed to run 
freely. Qij was set to zero throughout the grid every 10 timesteps for the first 
100 timesteps to take the simulation as close as possible to its minimal energy 
before the simulation proper began. 

3.6 Results 

The soliton-defect system displayed certain generic characteristics - a "spread- 
ing" channel and a "spiking" channel according to the phase of the soliton. In 
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the spreading channel the sohton would become broader and broader until it 
significantly overlapped the defect, at which point it would unwind, and the 
energy would be released to infinity as radiation. In the spiking channel the 
soliton would become more and more localised until the variation in the field 
became numerically untenable. The defect would remain stationary in all simu- 
lations, whilst the soliton would move only a very small distance, to the extent 
that variations in the width of the soliton would be far more significant than 
any movement of the maxima of the soliton. 

3.6.1 Stereographic Ansatz 

We used the stereographic ansatz (equation (|2^)) as our initial condition with 
A = 2.5, and the initial soliton position relative to the defect at 6d = tt, = 2.5. 
We shall discuss the results of simulations with 7 set to ^, — , and — j . After 
the initial period of relaxation energy was conserved to better than 0.05% in all 
of these simulations. The period of relaxation ended a,t t = 0.5, after which the 
simulation ran freely until t = 9.995, unless the simulation became untenable 
due to soliton spiking. 

When the simulation was run with the phase initially set to 7 = —-j the 
soliton became narrower with time as described above (see figure (^)). This 
resulted in the simulation ending at t = 3.745. 

When 7 = f the soliton spread out and eventually unwound as described 
above (see figure (^)). When 7 = was simulated the soliton was still subject 
to spreading and unwinding - the spreading channel is wider than the spiking 
channel. 

For 7 = — f the soliton began to spike, but after a period of time the 
the rate of spiking slowed and the soliton began to spread, with the soliton 
eventually unwinding (see figurc(^). This also suggests that the spiking channel 
is unstable. 

Table (|l|) shows the energies of various simulations with initial conditions 
given by this ansatz. Note that these energies depend heavily on the lattice spac- 
ing at the defect, and are therefore only of value when comparing the different 
channels. 
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Initial Energy 


Energy after Relaxation 


Time at end 


Final 




it = o) 


{t ^ 0.5) 


of Simulation 


Energy 


TT 

2 


1.439 


1.435 


9.995 


1.435 





1.520 


1.512 


9.995 


1.512 


TV 

2 


1.595 


1.590 


3.745 


1.591 



Table 1: Energy of System in Stereographic Ansatz 
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Figure 4: Energy density for 7 = | at t=0, 2, 4, 6, 8 and 9.975 
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Figure 5: Energy density for 7 = — ^ at t=0, 1, 2 and 3 
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Figure 6: Energy density for 7 = — f at t=0, 2, 4, 6, 8 and 9.75 
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3.6.2 Inhomogeneous Ansatz 



We also simulated the soliton defect system using the inhomogeneous ansatz 
(equation (|3l|)) as the initial condition, with the initial A = 2.5, and the initial 
soliton position relative to the defect at 9d — 0,rd — 2.5, so that the contribution 
to the ansatz from the soliton is significant compared to that from the defect in 
the region of the soliton. With this ansatz the soliton spikes for 7 = f (see figure 
(||)) and spreads at 7 = (see figure (^). If 7 = ^ in the initial condition, the 
soliton starts to spike, but then later spreads (see figure(||)). Table (||) shows 
the energies for simulations with these initial conditions. Again, these energies 
are only included for their value in comparing channels and to show that energy 
is conserved during the free run of the simulation. 
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Initial Energy 


Energy after Relaxation 


Time at end 


Final 




(t = 0) 


{t = 0.5) 


of Simulation 


Energy 


7r 
2 


1.597 


1.593 


3.745 


1.593 





1.523 


1.516 


9.995 


1.516 


TV 

2 


1.442 


1.437 


9.995 


1.438 



Table 2: Energy of System in Inhomogeneous Ansatz 



3.7 Collective Coordinate Approach 

To further the understanding of our numerical results we carried out some ap- 
proximate analytic work based on the so called collective coordinate approach. 
We substituted one of our ansatzc into the Lagrangian density, integrated over 
the spatial variables to find a true Lagrangian and then considered one or more 
of the ansatz parameters as dynamic variables, so for example Lagrangian (^ 
together with 

(f ^ (l){X,X,x,y), (33) 

implies that 

C = C{X,Kx,y), (34) 
which may then be integrated over to give 

L = y C{\, \,x,y)dxdy = L{\, X), (35) 

allowing us to construct an equation of motion for A. 

This effectively constrains the solutions of the field theory to the submanifold 
defined by the ansatz - in the example above our solutions are constrained to 
a 2 dimensional submanifold of the infinite dimensional phase space of (j){x,y). 
This does not allow for soliton unwinding or the release of radiation with our 
ansatze, as the ansatze never have radiation and always have a soliton. 
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Figure 7: Energy density for 7 = -f at t=0, 2, 4, 6, 8 and 9.975 
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Figure 8: Energy density for 7 = | at t=0, 1, 2 and 3 
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Figure 9: Energy density for 7 = -j at t=0, 2, 4, 6, 8 and 9.75 



19 



If we take the ansatz (|27|) together with the Lagrangian (|1|) expressed in 
terms of stereographic coordinates 

_ d^Wd>^W* 

4(1 + |W^P)2' W 

we may construct the Lagrangian density in terms of position and a few param- 



eters. If we take the initial sohton and defect position outhned in section (3.6.1) 
then W becomes 

_ 2A(rsin(6l-7)-2.5sin7) + i(l-A2(r2+5rcos6' + 6.25)) 9 
l + A2(r2 + 5rcos6' + 6.25) + 2A(rcos(6l-7) + 2.5cos7) ^^"^4' 

where (r, 6) are polar coordinates, A parameterises the (inverse) width of the 
soliton and 7 is the phase of the soliton. Note that the soliton has position 
(2.5, tt) and the defect is at the origin in this coordinate system. Taking spatial 
derivatives of W is then a straight forward if somewhat tedious process. Time 
derivatives may be found by treating one or more of the parameters as dynamic 
- if we consider A to be dynamic and the other parameters to be static the 
dt — \d\. This substitution of an ansatz and use of parameters as dynamic 
variables is equivalent to assuming that the field configuration moves quasi- 
statically from one configuration to another with different values of the dynamic 
variables. This assumption is valid only if our ansatz is close to equilibrium and 
our dynamic parameter only varies slowly with time (i.e. in this case A is small). 
Using this approximation our Lagrangian density becomes 

r^drWdrW* + deWdeW* - r^X^dxWdxW* 
^ 4r2(l + |iy|2)2 ■ ^ ' 

So our approximate Lagrangian with a time dependant A becomes 

A{\)->?B{X), (39) 

where 

^(^) - j 4r2(l + |Ty|2)2 ^'^^'^^ (40) 



and 



C d\Wd\W* 

B{\) = / "^"^'"^'' rdrde. (41) 

^ ^ J 4(1+|M^|2)2 ^ ^ 

The Euler-Lagrange equation then gives us 

■■ -iA'iX) + X^B'iX)) 

^' 2B{X) • ^^^> 

We considered the example above, where all variables are static except for 
A, but found an analytic integration to be intractable. We carried out the 
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integration and time evolution numerically, using a fourth order Runge-Kutta 
algorithm for the time evolution. We also carried out a similar analysis using 
the inhomogeneous ansatz (|3l|). 

In figures (jl^, |ll|) we show the results of this treatment against those of the 
full simulation. The two treatments produce broadly similar results, although 
the rapidity of the broadening and spiking is faster in the full simulation. This is 
not entirely unreasonable considering that we have moved from around 160,000 
degrees of freedom to 2! The curves from the full simulation were found by 
finding the maximum winding number density on the grid and then finding the 
value of lambda that would give a single soliton of the form of equation ( p4|) 
this maximum winding number density. One consequense of this technique is 
that for broad solitons the maximum winding number for the soliton may be 
smaller than the maximum on the grid - this leads to the curves for the broad- 
ening channel becoming unreliable at around t = 5. Naturally, the collective 
coordinate simulation began at the point when the relaxation ended in the full 
simulation. 




Figure 10: A vs. time for spiking and broadening channels in the stereographic 
ansatz 




Figure 11: A vs. time for spiking and broadening channels in the inhomogeneous 
ansatz 
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4 Conclusions 



We have examined the sigma model and baby Skyrme model with M.P'^ as a 
target space and found these models to be identical to their counterparts in 
the absence of defects. We examined the interaction between defects and soliton- 
like lumps in the sigma model and found the interaction to depend on the relative 
phase of the lump. We found a channel which causes the soliton to broaden 
and another which causes the soliton to spike. When the soliton overlapped 
the defect significantly, the soliton would unwind. We broadly reproduced this 
behaviour with the collective coordinate approach, using only two collective 
coordinates. 

A range of possibilities for future work present themselves - studying the 
interaction between two defects, the interactions of two lumps in the presence 
of a defect and the behaviour of the MP^ sigma model on a torus all have the 
potential to exhibit new and interesting behaviour. 
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